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Abstract Compared to the previous INPOP versions, the INPOPlOa planetary 
and lunar ephemeris has several improvements. For the planets of our solar system, 
no big change was brought in the dynamics but improvements were implemented 
in the fitting process, the data sets used in the fit and in the selection of fit- 
ted parameters. We report here the main characteristics of the planetary part of 
INPOPlOa like the fit of the product of the Solar mass with the gravitational con- 
stant (GM©) instead of the astronomical unit. Determinations of PPN parameters 
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as well as adjustments of the Sun J 2 and of asteroid masses are also presented. 
New advances of nodes and perihelia of planets were also estimated and are given 
here. As for INPOP08, INPOPlOa provides to the user, positions and velocities of 
the planets, the moon, the rotation angles of the Earth and the Moon as well as 
TT-TDB chebychev polynomials at http://www.imcce.fr/inpop. 

Keywords Planetary ephemerides • Celestial mechanics • Fundamental physics 

1 Introduction 

Since INPOP08 new observations have been produced, in particular normal points 
deduced from the three flybys of the Messenger spacecraft around Mercury. These 
data were the first accurate positions of Mercury since the two Mariner flybys 
in the seventies. Furthermore, the asteroid masses obtained during the fit of IN- 
POP08 (Fienga et al. 2009 [11 ) to observations did not seem to be satisfactory. 
The method used for these estimations has been revised and a more constraining 
approach was tested. The question of fitting the mass of the Sun instead of the 
astronomical unit in the planetary ephemerides was discussed in the community. 
We decided to test the possibility of such adjustment by building an ephemerides 
with a fixed value of AU and a fitted value of GM©. Finally, crucial normal points 
for the outer planet orbits missed in the INPOP08 data sets have been added. 

In this paper, we introduce the planetary ephemerides INPOPlOa and some 
applications. As no modification was brought to the dynamical modeling, we first 
describe the main input to the observational data sets used for the adjustment. We 
then introduce the modifications implemented in the selection of used constants 
and fitted parameters. In the section 4, we describe the obtained results in term 
of comparisons to other planetary ephemerides (INPOP08, DE421 (Folkner et 
al. 2008 [H])) and of postfit comparisons to observations. The first combined 
estimations of PPN parameters f3 and 7 based on INPOP are presented here as 
well as a list of 58 asteroid masses estimated with INPOPlOa. Supplementary 
advances of perihelia and nodes of planets are also given. Finally we introduce 
the first estimation of rotation angles between planetary ephemerides frames and 
the ICRF deduced from the VLBI and radio astrometry of millisecond pulsars. We 
check the internal mas-level accuracy of the Earth orbit but we detect a possible 10 
milliarcseconds (mas) rotation between all the planetary solutions and the ICRF. 
This result needs to be verified after the densification of pulsars with a mas-level 
VLBI and radio astrometry. 

2 The INPOPlOa Data Sets 

A detailed description of the observations used for the construction of INPOPlOa 
can be found in (Fienga et al. 2010 [H]). Several data sets have been added since 
INPOP08. The global observational distribution for the INPOP fit has changed its 
balance compared to INPOP06: now, more than 56% of the planetary observations 
are deduced from the tracking data of spacecrafts including range, VLBI angular 
positions and flyby normal points. The statistics of the obtained postfit residuals 
as well as the number of points and their distribution in time are presented in 
table □ 
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The first spacecraft flying by Mercury was Mariner 10. From this mission and 
its two flybys of Mercury in 1974 and 1975 were extracted two corrections to 
the geocentric distance of the planet, provided by JPL (Folkner 2010 [14 ) with 
an accuracy of about 150 meters. The MErcury Surface, Space ENvironment, 
GEochemistry, and Ranging (MESSENGER) spacecraft, only the second probe to 
encounter the innermost planet, flew by Mercury on 14 January 2008, 6 October 
2008, and in December 2009. 

The MESSENGER payload includes the Mercury Laser Altimeter (MLA) (Ca- 
vanaugh et al. 2007 [5]) and a gravity science investigation that uses the spacecraft 
tracking system (Srinivasan et al. 2007 [38]). During the flybys the spacecraft was 
tracked at X-band frequency by the Deep Space Network (DSN), the latter provid- 
ing high-quality Doppler range-rate observations of the spacecraft motion relative 
to Earth. From the three flybys and the accurate orbit determination obtained 
by (Smith et al. 2010 [37 J, we extracted three corrections to Mercury geocen- 
tric positions with an accuracy of about a few mas in angles (right ascension and 
declination) and a few meter in Earth-Mars range. The MESSENGER orbit recon- 
structed from tracking data by the mission navigation team is available through the 
NASA Planetary Data System and can be accessed with SPICE software (Turner 

2007 06]). With the Mercury-centric orbit of the spacecraft, the navigation team 
provides also orbits of planets (Earth, Mercury and Venus) fitted only during the 
flybys of the spacecraft. We then deduced for each flyby of Mercury, corrections 
to Mercury's geocentric orbit in angular positions and distances. 

These five points changed our knowledge of the Mercury orbit: until now, only 
direct radar ranging on Mercury's surface were available with an accuracy of about 
800 meters. 

For Mars and Venus, like with INPOP08, tracking data of MEX and VEX 
missions provided by ESA (Morley 2009 [29], Morley 2010 [30 J are used in the fit 
as described in Fienga et al. (2009) [11] , In addition to the Saturn Cassini normal 
points provided by JPL over the 2005 to 2007 period and used in INPOP08, are 
added VLBI observations of the spacecraft (Jones et al. 2011 [18]) with an accuracy 
better than a few mas. These VLBI differential positions of Cassini related to ICRF 
sources were obtained during its mission about Saturn and its satellites. 

Flybys data of Jupiter, Uranus and Neptune obtained during several missions 
(Pioneer 10 and 11, Viking 1 and 2, Ulysses and Cassini) were also added, provided 
by Folkner (2009) [13]. These observations improve the estimates of the geocentric 
distances of the outer planets while no observations of similar type were used in 
INPOP06 and INPOP08 adjustments. New optical data obtained from 2000 to 

2008 with the Flagstaff Astrometric Scanning Transit Telescope are also added 
for Uranus, Neptune and Pluto (Folkner 2009 [13]). Stellar occultations (Sicardy 

2009 [36]) are taken into account in INPOPlOa by the use of measured offsets in 
topocentric (a,S). 



3 INPOPlOa General features 

In INPOPlOa, the motion of the planets including Pluto and the Moon are inte- 
grated with the rotations of the Earth and the Moon, the (TT-TDB) differences 
and the orbits of several hundreds asteroids as in INPOP08 (Fienga et al. 2009 
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Table 1 Statistics of the residuals obtained after the INPOPlOa fit. For comparison, means 
and standard deviations of residuals obtained with INPOP08 are given. 



Planet 






INPOP08 


INPOPlOa 


Type of Data 


Nbr 


Time Interval 


mean 


la 


mean 


la 


Mercury- 














Direct range [m] 


462 


1965-2000 


30 


842 


7 


866 


Mariner range [m] 


2 


1974-1975 


-1000 


305 


-28 


85 


Messenger Mercury flybys 














Messenger ra [mas] 


3 


2008-2009 


1.1 


0.7 


0.4 


1.2 


Messenger de [mas] 


3 


2008-2009 


2.0 


1.9 


1.9 


2.1 


Messenger range [m] 


3 


2008-2009 


52 


619 


-0.6 


1.9 


Venus 














Direct range [km] 


489 


1965-2000 


0.5 


2.3 


0.5 


2.2 


VEX range [m] 


22145 


2006-2010 


1.6 


4.4 


-0.2 


3.9 


VLBI [mas] 


22 


1990-2007 


2 


2 


2 


2.5 


Mars range [m] 














MGS 


10474 


1998-2008 


-0.9 


1.6 


0.5 


1.9 


MEX 


24262 


2006-2010 


-3.5 


2.0 


0.0 


1.7 


Path 


90 


1997 


6.8 


12.5 


-5.0 


5.0 


Vkg 


1256 


1976-1982 


-27.4 


19.0 


-5.7 


35.0 


Mars VLBI [mas] 


96 


1989-2007 


0.5 


0.5 


-0.0 


0.4 


Jupiter flybys 














ra [mas] 


5 


1974-2000 


48.0 


40.0 


6 


5 


de [mas] 


5 


1974-2000 


-10.0 


50 


-13 


18 


range [km] 


5 


1974-2000 


-27 


55 


-0.6 


1.6 


Jupiter VLBI [mas] 


24 


1996-1997 


4 


11 


0.2 


11 


Jupiter Optical 














ra [mas] 


6216 


1914-2008 


20 


304 


-26 


304 


de [mas] 


6082 


1914-2008 


-44 


313 


-54 


303 


Saturn Cassini 














ra [mas] 


31 


2004-2007 


1.5 


4 


0.7 


4 


de [mas] 


31 


2004-2007 


7.0 


7 


6.5 


7 


range [m] 


31 


2004-2007 


0.5 


22 


0.0 


17 


Saturn VLBI Cassini 














ra [mas] 


10 


2004-2009 


0.3 


0.7 


0.0 


0.6 


de [mas] 


10 


2004-2009 


-1.2 


2.0 


0.1 


0.4 


Saturn Optical 














ra [mas] 


7824 


1914-2008 


-16 


305 


-16 


305 


de [mas] 


7799 


1914-2008 


-7 


276 


-9 


276 


Uranus flybys 














ra [mas] 


1 


1986 


-90 




-30 




de [mas] 


1 


1986 


-36 




-7 




range [km] 


1 




1139 




0.080 




Uranus Optical 














ra [mas] 


4145 


1914-2008 


-44 


278 


-27 


290 


de [mas] 


4130 


1914-2008 


-38 


339 


-11 


338 


Neptune flybys 














ra [mas] 


1 


1989 


-88 




-11 




de [mas] 


1 


1989 


-48 




-10 




range [km] 


1 




2305 




0.004 




Neptune Optical 














ra [mas] 


4340 


1914-2008 


-32 


282 


2 


281 


de [mas] 


4320 


1914-2008 


-36 


335 


2 


330 


Pluto occultation 














ra [mas] 


13 


2005-2009 


-6 


46 


-1 


47 


de [mas] 


13 


2005-2009 


16 


30 


-2 


19 


Pluto Optical 














ra [mas] 


2449 


1914-2008 


353 


926 


38 


629 


de [mas] 


2463 


1914-2008 


-22 


524 


17 


536 
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Table 2 Asteroid masses found in the recent literature and compared to the values estimated 
in INPOPlOa. [1]: (Baer et al 2011 0), [2]: (Baer et al 2008 OQ), [3]: (Marchis et al. 2008a [26]), 
[4]: (Marchis et al. 2008b [27]), [5] : (Yeomans et al. 1997 48 )", [6]: (Kochetova 2004 [19]). Kll 
stands for (Konopliv et al. 2011) 20 . [*]: INPOP08e. The last column gives the impact of each 
asteroid on the Earth-Mars distances over the 1970 to 2010 period. The impact is defined as the 
difference in meters between two integrations of the Mars geocentric orbit, with and without 
the corresponding asteroid. The uncertainties for INPOPlOa and INPOP08 are given at 1 
formal sigma deduced from least squares. For the close-encounters and Kll determinations, 
the uncertainties are given at 1 published sigma 



Asteroids INPOPlOa Others INPOP08 Kll Impact 

1970-2010 



IAU 


1 1 2 , 71 /T 

10 X Mq 


1 /~1 1 2 ^ , 7\/T 

10 X Mq 




1 n1 2 , 7\/T 
10 X Mq 


10 X Mq 


m 


1 


475.8 zb 2.8 


475.700 zb 0.72 


[!] 


465.8 zb 0.9 


467.900 zb 3.250 


794 


2 


111.4 ± 2.8 


101.000 ± 6.5 


[1] 


107.6 ± 2.0 


103.440 ± 2.550 


146 


4 


133.1 zb 1.7 


130.00 zb 0.53 


[1] 


132.9 zb 3.0 


130.970 zb 2.060 


1199 


7 


7.7 zb 1.1 


8.12 zb 0.46 


[1] 


5.0 zb 0.2 


5.530 zb 1.320 


28 


324 


4.67 zb 0.38 






5.6 zb 0.08 


5.340 zb 0.990 


94 


3 


11.6 zb 1.3 


14.400 zb 2.3 


[1] 


7.5 zb 0.3 


12.100 zb 0.910 


56 


6 


7.1 zb 1.2 


6.40 ± 0.67 


[1] 


0.16 ± 0.11 


6.730 ± 1.640 


21 


8 


4.07 zb 0.63 


3.33 ± 0.42 


[1] 


5.35 zb 0.5 


2.010 zb 0.420 


13 


9 


5.700 [1] 


5.700 ± 1.1 


[1] 


1.15 zb 1.0 


3.280 zb 1.080 


30 


10 


44.500 [2] 


43.58 ± 0.74 


[1] 




44.970 zb 7.760 


77 


11 


1.9 zb 1.0 


3.090 ± 0.989 


[2] 






17 


13 


8.200 [2] 


8.00 ± 2.2 


[1] 






<1 


14 


4.130 [2] 


3.49 ± 0.82 


[1] 


0.71 zb 0.1 


1.910 zb 0.810 


18 


15 


18.8 zb 1.6 


15.597 ± 0.15 


[1] 


22.30 ± 1.80 


14.180 ± 1.490 


22 


16 


11.2 zb 5.2 


11.40 ± 0.42 


[1] 


15.96 zb 0.32 


12.410 zb 3.440 


10 


18 


1.845 [2] 








11 


19 


6.380 [2] 


4.18 ± 0.36 


[1] 


2.02 zb 0.20 


3.200 zb 0.530 


59 


20 


2.850 [2] 


1.68 ± 0.35 


[1] 






< 1 


21 


1.3 zb 1.2 


1.31 ± 0.44 


[1] 


1.034 zb 0.30 




5 


24 


2.8 zb 1.9 


5.670 ± 2.155 


[2] 






26 


25 


0.002 zb 0.002 






0.301 zb 0.030 




3 


28 


4.65 zb 1.0 










4 


29 


5.920 [2] 


7.63 ± 0.31 


[1] 


4.91 zb 0.90 


7.420 zb 1.490 


27 


31 


3.130 [2] 


2.92 ± 0.99 


[1] 


29.9 zb 6.8 




23 


41 


9.2 zb 2.6 






5.27 zb 0.5 


4.240 zb 1.770 


12 


42 


1.8 zb 1.07 










7 


45 


2.860 [3J 


2.860 ± 0.057 


To! 

[3] 






<1 


51 


0.72 zb 0.42 










5 


52 


42.3 ± 8.0 


11.39 ± 0.79 


[1] 


17.24 ± 1.0 


11.170 ± 8.400 


10 


60 


0.402 ± 0.37 










6 


63 


2.022 zb 1.68 










6 


65 


7.2 zb 4.2 


p*" o n i r\ f\r* 

5.30 ± 0.96 


[i] 






5 


78 


3.23 zb 2.3 










9 


94 


15.8 ± 11.5 










7 


105 


1.109 [*] 










<1 


1 07 
1U / 


lo.Z ZTZ 4t.D 


o.uou zc u.ioy 








K 



130 


11.1 ± 8.0 


3.320 ± 0.199 


[4] 






<1 


135 


0.917 ± 0.88 










2 


139 


5.9 zb 3.3 






3.59 zb 0.10 




17 


145 


2.266 [*] 










<1 


187 


2.5 zb 1.07 










<1 


192 


0.719 [*] 






1.37 zb 0.41 




11 


194 


8.8 zb 2.9 










5 


216 


0.560zb 0.46 






3.53 zb 0.12 




3 


253 


0.904 zb 0.65 


0.052 zb 0.002 


[5] 






<1 


337 


0.543 ± 0.080 










2 


344 


0.340 zb 0.19 










7 


354 


2.451 [*] 






4.88 zb 0.35 




10 


372 


4.443 [*] 










<1 


419 


0.997 zb 0.55 










10 


451 


21.0 ± 14.8 


10.2 zb 3.4 


[6] 






5 


488 


6.2 zb 5.5 










9 


511 


19.9 ± 4.1 


18.96 zb 0.99 


[1] 




8.580 zb 5.930 


10 


532 


2.89 zb 0.76 


16.8 zb 2.8 


[6] 


5.46 zb 0.10 


4.970 zb 2.810 


34 


554 


1.6 zb 1.3 










5 


704 


18.600 [2] 


19.65 zb 0.89 


[1] 


16.23 zb 1.0 


19.970 zb 6.570 


34 


747 


6.0 zb 2.3 






0.04 zb 0.02 




16 


804 


2.5 zb 1.8 


1.75 zb 0.40 


[1] 
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Table 3 Values of parameters obtained in the fit of INPOP08 and INPOPlOa to observations. 
The (F) indicates that the marked values are fixed in the fit. The equivalent value of AU 
deduced from the estimation of the GM© in INPOPlOa and (Konopliv et al. 2011) [20| are 
given in the line labelled "AU from GM0" . For the AU, the GM© and the asteroid masses, 
the uncertainties are deduced from the least squares fit. 





INPOP08 
± la 


INPOPlOa 
± la 


Kll 


(EMRAT-81.3000)x 10~ 4 


(5.4 ± 0.5) 


(5.7 ± 0.010) 


(5.694 ± 0.015) 


J 2 ° x 10- 7 


(1.82 ± 0.47) 


(2.40 ± 0.25) 




GM - GMq E405 [km 3 , s- 2 ] 


0.0 ± 50 (F) 


37.013 ± 1 


24.013 ± 10 


AU-AU /MS03 [m] 


8.2 db 0.11 


0.0 (F) 




AU [m] from GM© 




13.9 ± 0.3 


9.0 ± 3 



Table 4 Values of parameters obtained in the fit of INPOP08 and INPOPlOa to observa- 
tions. The given values are (3 and 7 intervals in which the differences of postfit residuals from 
INPOPlOa are below 5%. Kll stands for (Konopliv et al. 2011) 20 , M08 for (Muller et al. 
2008) [31], W09 for (Williams et al. 2009) 47 , B03 for (Bertotti et al. 2003) 4 and L09 for 
(Lambert and al. 2009) | 23] 

Fixed PPN Estimated PPN x 10 4 



7 


= 1 


(P 


-1) 


INPOP08 


(0.75 ± 1.25) 










INPOPlOa 


(-0.62 ± 0.81) 


7 


= 1B03 


(P 


-1) 


INPOPlOa 


(-0.41 ± 0.78) 










Kll 


(0.4 ± 2.4) 










M08-LLR-SEP 


(1.5 ± 1.8) 










W09-LLR-SEP 


(1.2 ± 1.1) 


P 


= 1 


(7 


-1) 


INPOPlOa 


(0.45 ± 0.75) 










Kll 


(1.8 ± 2.6) 










B03-CASS 


(0.21 ± 0.23) 










L09-VLB 


(0.7 ± 2.0 ) 



[TT]). The initial conditions of the planet orbits are fitted over the data sets de- 
scribed in section [5] The Moon orbit and rotation initial conditions are fitted 
separately by Manche et al. (2010) [25]. An iterative process between planetary 
and lunar adjustments is used to keep consistencies between the two part of the 
ephemerides. With the initial conditions of the planets are also fitted the Earth- 
Moon mass ratio and the oblateness of the Sun coefficient J2. 

In opposition to INPOP08 where the value of astronomical unit (AU) is fitted, 
in INPOPlOa, the value of the AU is fixed to the value given by the IERS2003 
conventions [28]. The GM of the Sun is fitted to the observations like the other 
fitted quantities. Values of the fitted parameters are given in table [3] (Konopliv 
et al. 2011) [20] presents also a value of the GM of the Sun. The two estimations 
are consistent at 2 sigmas with a ten times smaller uncertainty for the INPOPlOa 
estimation. This can be explained by the fact that the INPOPlOa uncertainties 
are 1 formal sigma deduced from least square fit when the uncertainties given by 
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Table 5 Values of parameters obtained in the fit of INPOP08 and INPOPlOa to observations. 
The supplementary advances of perihelia and nodes are estimated in INPOPlOa and INPOP08 
as the interval in which the differences of postfit residuals from INPOPlOa are below 5%. P09 
stands for (Pitjeva 2010 [34]). 





TNPOPOS 


TNPDP1 Da 


P09 










mas x cy — 1 








Mercury 


-1U zb oU 


U.4 zb U.O 


-4 zb 


^^enus 


-4 zb 6 


0.2 zb 1.5 


24 zb 33 


EMB 


0.0 zb 0.2 


-0.2 zb 0.9 


6 zb 7 


Mars 


0.4 zb 0.6 


-0.04 zb 0.15 


-7 zb 7 


Jupiter 


142 zb 156 


-41 zb 42 


67 zb 93 


Saturn 


-10 zb 8 


0.15zb 0.65 


-10 zb 15 


^sup 








mas x cy — 1 








Mercury 




1.4 zb 1.8 




Venus 


200 zb 100 


0.2 zb 1.5 




EMB 


0.0 zb 10.0 


0.0 zb 0.9 




Mars 


0.0 zb 2 


-0.05 zb 0.13 




Jupiter 


-200 zb 100 


-40 zb 42 




Saturn 


-200 zb 100 


-0.1 zb 0.4 





(Konopliv et al. 2011) [20] also involved the propagation of the error induced by 
the use of fixed asteroid masses in the fit. 



3.1 The fit procedure for the asteroid masses 

The estimations of the asteroid perturbations on planet orbits are a critical point 
for the extrapolation capabilities of the planetary ephemerides (Standish and 
Fienga 2002 [39]). 

The usual approach to this problem has been suggested by Williams 1984 and 
consists of including for a selection of approximately 300 individual asteroids to 
the dynamical model. The masses of the most perturbing asteroids are fitted to 
observations. For the other objects, masses are deduced from radiometric diameters 
and the assumption of constant densities within three taxonomic classes. This 
classic approach has been used in INPOP08 and achieves in terms of the Earth- 
Mars distance prediction an accuracy of 20 m over 2 years. It is based on an 
unrealistic hypothesis of constant densities within taxonomic classes. It also relies 
on an empirical choice of the selection of asteroids to account for and on the choice 
of the subset of asteroid masses to adjust individually. We used INPOPlOa as a 
benchmark to test an alternative approach. Kuchynka et al. (2010) [22] showed 
that approximately 240 asteroids in a list of 287 probable asteroids and a ring 
should represent the perturbations induced by the main belt on planetary orbits 
down to an order of a meter. The result was obtained for various test models 
of the main belt containing each several thousands of asteroids with randomly 
assigned masses. In INPOP10, we use the Bounded Variable Least Squares (BVLS) 
algorithm developed by (Lawson and Hanson 1995 [24]) and (Stark and Parker 
1995 02]) in order to fit the masses of all the 287 asteroids listed in Kuchynka et 
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al. (2010) [22]. The BVLS solves the least square problem with constraints which 
require the adjusted masses to be positive or zero. Setting an asteroid mass to zero 
is equivalent to removing it from the dynamical model. Thus the BVLS algorithm 
performs simultaneously parameter selection and estimation. From this method 
and the original list of 287 asteroids, about 161 asteroid masses were planed to be 
estimated in the fit, the other masses being put to zero. 

Based on a study of the correlations between the asteroid masses, we found 
however 30 highly correlated asteroids among the most perturbing objects. In 
order to decrease the mass uncertainties, we fixed 16 asteroid masses to values 
determined by other methods (close encounters, binary) or by previous published 
or unpublished INPOP versions as INPOP08 or INPOP08e. in table 2, the fixed 
masses are identified with their references. 

Besides these fixed values, we then estimate 145 asteroid masses using con- 
straints on the densities in order to keep the fitted values in the frame of re- 
alistic physics. These densities are constrained to be positive and smaller than 
20 gxcm -3 . The diameters used to estimate the densities are extracted from a 
database (Kuchynka 2010) [21 compiling mainly diameters from the IRAS cata- 
logue (Tedesco et al. 2002 [45]). For objects without radiometric data, diameters 
are estimated from assumed mean albedo (Tedesco et al. 2005 04]). 



4 Results and applications 

4.1 Postfit residuals and comparisons with other planetary ephemerides 

Complete plots of comparisons between INPOP06, INPOP08, INPOPlOa and 
DE421 can be find in Online Resource 1 (Figures 1 to 6). For all the planets, 
one can note the convergence of the last INPOP solution and DE421. We stress in 
particular the reduction of the differences between planetary ephemerides for the 
Jupiter orbit. As one can see in figure 4, INPOP08 appears then to be biased es- 
pecially for Jupiter orbit. This can be explained by the fact that no normal points 
deduced from mission flybys were included in INPOP08 while a large amount of 
Mars mission data were presented. The INPOP08 estimated orbit of Jupiter seems 
then to be over-weighted by Mars data. In the case of INPOP06, no flyby points 
were included but the fit was far less Mars observation-depend than INPOP08. 
In INPOPlOa, Mars observations are still very important but the input of the 
flyby normal points balances the impact for Jupiter's orbit. Furthermore, due to 
transmission problems during the mission, VLBI Galileo angular positions have 
a degraded accuracy compared to the one expected . Thus, even if these points 
were included in the INPOP08 adjustment, they were not reliable enough to cor- 
rect Jupiter's orbit. For Venus, the input of the VEX data is obvious in figure 1 
as for Saturn and the impact of Cassini data. The noticeable offsets for Saturn 
and Jupiter longitudes between INPOPlOa and DE421 are far below the limit of 
accuracy of the data used for the fit. For Uranus and Neptune (figure 5), the dif- 
ferences between the ephemerides are about the same with a change of behavior 
induced by the input of the flyby points in 1986 for Uranus and 1989 for Neptune. 
For Mercury, despite the input of the five flyby points, the differences between the 
ephemerides including or not these points are about the same. 
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Table [T] gives the planetary postfit residuals obtained with INPOP08 and IN- 
POPlOa. Some of the data sets used in INPOPlOa were not used in the INPOP08 
adjustment: this explains some of the differences between the two columns of Ta- 
ble [l] Figures 1 to 6 in Online Resource 2 and Figures 1 to 3 in Online Resource 
3 present the postfit residuals obtained with INPOPlOa for the different sets of 
observations used in the fit. As one can see with the extrapolated residuals of 
INPOP08 for the outer planet flybys (see Table [TJ, the improvement of the inner 
planets came together with a degradation of the outer planet orbits, no noticeable 
with the postfit residuals obtained at that time (Fienga et al. 2009 [TT]). In other 
words, the outer planet orbits of INPOP08 are accurate at the level of uncertainty 
of the optical and VLBI observations. INPOPlOa, with the use of the flyby data 
of the outer planets and the improvement of the asteroid masses, provides good 
residuals for all the planets. Figures of INPOP06 and INPOP08 residuals obtained 
with the INPOPlOa data set are given in Online Resources 2 and 3. 



4.2 Asteroid masses 

in figure [l] are presented the asteroid densities deduced from INPOPlOa compared 
to values found in the literature, ranked by their impact on the Earth-Mars dis- 
tances over the 1970 to 2010 period. The major sources of comparisons are the 
values obtained with INPOP08, (Konopliv et al. 2011) [20] and close encounters 
or binary system estimations gathered by (Baer et al. 2011 [2 ). The densities are 
deduced from the mass estimations and diameters extracted from the (Kuchynka 
2010) database. 

It appears clearly from figure [l] that the estimations for the most perturbing 
objects are quite consistent while the estimations of the less perturbing objects 
show bigger discrepancies. 

in tables[2]are given the details of asteroid masses obtained with the adjustment 
of INPOPlOa. As expected the three bigger asteroid masses are compatible at 2- 
sigmas. 

in figure [2] are plotted the distributions of densities deduced from INPOPlOa, 
INPOP08 and from close encounters and binary systems (noted "Other"). Two 
representations are given: one histogram of density distribution and one distri- 
bution of the density versus the diameters of the objects. Are plotted on these 
figures, only the densities deduced from perturbations bigger than 1 meter on the 
Mars-Earth distances over the 1970 to 2010 period with error bars representing 
the 1-sigma uncertainties on the mass determination. The diameters are consid- 
ered here as perfect. With this optimistic hypothesis, one can first note the smaller 
uncertainties on the close-encounter estimations compared to those obtained with 
INPOP08 and INPOPlOa. The distribution of the close-encounter densities shows 
also two very close peaks of frequency at 2.5 and 3.5 g x cm -3 , and about few per- 
centages of objects with a density below 1 g x cm -3 . This trend is opposite to the 
INPOP08 distribution with more than 20 % of low density asteroids. This behavior 
was noticed in (Fienga et al. 2009 [TT]) and shows the excess of under-estimated 
masses in INPOP08. In the other hand, the dispersion and the uncertainties of the 
INPOPlOa distribution do not allow to give conclusive remarks even if one can 
note a diminution of the number of low density objects compared to INPOP08. 
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Fig. 1 Distribution of densities in g x cm -3 ranked with the impact on the Mars- Earth dis- 
tances over 1970-2010 period. The right-hand side axis gives the differences in Mars geocentric 
distances in meters induced by an integration of the Mars motion over the 1970-2010 period 
with and without the corresponding asteroid. Kll stands for (Konopliv et al. 2011) [20] and 
Other for the determinations induced by close-encounters or binary systems. The error bars 
represent the 1-sigma error on the mass determinations. 



4.3 Gravity tests 

Table [3] gives the obtained values for the GM of the Sun, the Sun J2, the Earth- 
moon mass ratio and table [2] the interval of sensitivity of data to modifications 
of PPN ft with 7 equal to 1 or with 7 equal to the value obtained by the Cassini 
experiment (Bertotti et al. 2003 |4 ), and of PPN 7 with f3 equal to 1. The results 
obtained for both modified parameters are presented in figure [3] Supplementary 
advances in the perihelia and nodes of the planets (from Mercury to Saturn) have 
also been tested and are presented in table [5] All these results were based on the 
method presented in Fienga et al. (2010) [10]. We computed several fits for different 
values of the PPN parameters (/3, 7) or supplementary advances in perihelia and 
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Fig. 2 Histograms of densities in g x cm -3 for the most perturbing objects deduced from 
INPOPlOa, INPOP08, and Other determinations (close encounters and binary system) and 
distributions of densities versus the diameters in kilometers. The plotted uncertainties represent 
the 1-sigma error on the mass determinations. 



nodes with a simultaneous fit of initial conditions of planets, mass of the Sun 
and asteroid densities. The given intervals correspond to values of parameters 
inducing changes in the postfit residuals below 5% compared to INPOPlOa postfit 
residuals, in figure [3] are plotted the variations of postfit residuals induced by 
the modification of the corresponding /3 and 7 parameters. The left hand side 
plot gives the variations of postfit residuals including Mercury flyby normal point 
when the right hand side plot gives the variations of residuals without the Mercury 
observations. 

The different levels of grey indicate the percentage of variations of the postfit 
residuals compared to those obtained with INPOPlOa. By projecting the 5% area 
on the /3-axis (or the 7-axis), one can deduced the corresponding f3 (or 7) interval 
in which the residuals are modified by less than 5%. In looking at the two figures, 
one can see that the use of the Mercury flyby data give smallest intervals of possible 
/3,7- This is consistent with the fact that the Mercury observations are far more 
sensitive to gravity modifications than other data (see table 1 in (Fienga et al. 
2010) pig ). 

As one can see in tableland in figure [3] the possible variations acceptable for 
all the residuals including the Mercury points are compatible with the estimations 
of /3 deduced from planetary ephemerides as (Konopliv et al. 2011) [20] and (Pitjeva 
2009) [34 . However, one can see that the /3 determinations constrained by Mercury 
flybys are not consistent with the one deduced from LLR measurements of the 
equivalence principle ((Muller et al. 2008) [31] and (Williams et al. 2009) [47]). If 
the Mercury flyby data are not taken into account (figure [I] right hand side plot), 
the constraint is then released and all the determinations become then compatible. 
The Messenger spacecraft is now orbiting Mercury. The use of its tracking data 
will help to improve the f3, 7 estimations and to confirm or not these denoted 
discrepancies. For the advances of perihelia and the nodes, the estimations based on 
INPOPlOa show the clear incompatibility of a significant supplementary advance 
in the planets' orbits and the observations used in the INPOPlOa adjustment. 
This conclusion is due to the densification of very accurate Cassini observations 
around Saturn. 
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Fig. 3 Variations of postfit residuals obtained for different values of PPN (5 (x-axis) and 7 (y- 
axis). [1] stands for a PPN (3 value obtained by (Manche et al 2010 [25]) using LLR observations 
with 7 = 0, [2] stands for Pitjeva (2009) 34 by a global fit of EPM planetary ephemerides. 
Kll stands for (Konopliv et al. 2011) 20 determinations based mainly on Mars data analysis. 
M08 for (Muller et al. 2008) [31] and W09 for (Williams et al. 2009) [47] give values deduced 
from LLR for a fixed value of 7, B03 stands for (Bertotti et al. 2003) 4 determination of 7 
by solar conjunction during the Cassini mission. 



4.4 Use of millisecond pulsars for reference frame ties 

One possible method to test the real accuracy of planetary ephemeris is to use very 
accurate observations of objects but not used in the fit process of the ephemerides. 
Nowadays three types of observations have a mas level astrometry: the VLBI 
surveys of ICRF sources, the VLBI tracking of spacecrafts into ICRF used in 
the construction of the planetary ephemerides and the millisecond pulsar (MSP) 
observations by radio timing and by VLBI astrometry relative to ICRF sources. 
Pulsar timing is the regular monitoring of the rotation of the neutron star by 
tracking (nearly exactly) the times of arrival of the radio pulses. By their relatively 
stable rotations and sharp pulses of nearly point-source radio emission, millisecond 
pulsars can have the astrometry accurate at the mas level if series of observations 
are done on a long enough time interval. 

Every single rotation of the neutron star is unambiguously accounted over long 
period of time. This unambiguous and very precise tracking of rotational phase 
allows pulsar astronomers to probe the physics interior of neutron stars, to make 
extremely accurate astrometric measurements, and test gravitational theories in 
the strong-field regime in unique ways. Of course, systematics in the pulse ob- 
servations exist such as the fluctuations of the interstellar medium (DM) or the 
binary structure of the neutron star. However, by having a long timespan of ob- 
servations it becomes possible to model the frequencies induced by the companion 
or to average the DM effect and then reduce significantly the error. Multi wave- 
length observations of millisecond pulsars are also performed in order to study the 
physics and the dynamics of these objects as the VLA polarization survey done by 
(Taylor et al. 1993 [43]) and (Han and Tian 1999 [17 ) or the several year program 
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using VLB A and VLBI techniques begun by (Chatterjee et al 2001 1 7 ) to estimate 
proper motions and parallaxes of about 100 pulsars. The mas level astrometry of 
pulsars are then an interesting tool for testing planetary ephemerides. By the com- 
bination of millisecond pulsar timing and VLBI observations of the same objects, 
it is possible to establish a new independent link between ICRF and dynamical ref- 
erence frames. (Shapiro and Knight 1970 [35]) have the first proposed to establish 
such direct link. 

Since 2009, VLBI observations of MSP are regularly done by (Chatterjee et 
al. 2009 [6]) in the northern hemisphere and (Deller et al. 2009 [8 ) in southern 
hemisphere. These two surveys provide positions of MSP in the ICRF at the mas 
level and give then a good opportunity to realize the link between ICRF and 
planetary ephemerides reference frames. 

VLBI-derived positions of MSP can provide a tie between the extragalactic ( 
Earth-rotation based) reference frame in which VLBI operates and the dynamic ( 
Earth-orbit based) reference frame in which MSP timing positions are derived. 

If {ot-PEifipEi) are positions of a MSP deduced from radio timing in using 
Planetary Ephemerides 1 (PE1) and {olpe2^pe2) are positions of the same MSP 
deduced either from VLBI observations in ICRF either from radio timing in using 
Planetary Ephemerides 2, the differences (ctpEifipEi) and (apE2,5pE2) can be 
seen as residual rotations R x (6),R y (rf) and R z (() about the x,y and z axis of PE1 
reference frame such as 

apEi = R x (0)R y (r])R z (()apE2 (1) 

and 

Spei = Rx(0)R y (ri)R z (()5pE2 (2) 

. The angles #,77 and £ can then be deduced by least square fitting. 

From the radio timing profiles obtained at the NRT (Desvignes 2009 [9 ), we 
selected 18 MSP with radio timing astrometry better than 10 mas. These 18 MSP 
presented in table [6] are used, as a check of the procedure, for the estimation of 
the rotation matrices between DE421 (Folkner et al. 2008 [16 ), INPOP08 (Fienga 
et al. 2009 [11]), DE405 (Standish 1998 [41]) and DE200 (Standish 1990 00]) 
presented in table [7] For the second step, we collected in the literature positions 
and velocities of millisecond pulsars obtained by VLBI astrometry (Chatterjee et 
al. 2009 ©, Deller et al. 2009 [8]) and observed by the NRT. 4 MSP have a mas- 
level accuracy in both techniques, VLBI and radio timing, and are used to test the 
link between planetary ephemerides frames and ICRF. Results are presented in 
table [7] 

We have first estimated the impact of using different planetary ephemerides 
(with different sets of planet masses) on the analysis of radio timing data in order 
to check the method. For the 18 MSP used for this study, no specific trend remains 
in the postfit residuals of the radio timing despite the change in planet masses and 
initial conditions brought by each different ephemerides. Among the parameters 
fitted during the radio timing data analysis, the positions and the proper motions 
of the pulsars are sightly modified. No other parameters related to the distance, 
the rotation or the orbit of the companion are affected by the change of planetary 
ephemerides. 

Based only on radio timing observations, the angles presented in table [7] and 
obtained in using the 18 MSP are quite compatible with previous determinations 
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Table 6 Main characteristics of millisecond pulsars used for this study. The last column 
indicates by a x if the pulsar is used for the ICRF link. [6] stands for (Chatterjee et al. 2009), 
8 for (Deller et al. 2009), [32] for (Nunes and Bartel 1998), [3] for (Bartel et al. 2006), [33] 
for (Petit 1994). 







Radio Timing 


VLBI 




MSP 


RMS 


cr a 




Period 


Ref 






/is 


mas 


mas 






mas 


J0139+58 


62 


800 


70 


2005-2009.5 


M 


15 


J0454+55 


90 


100 


100 


2004.5-2009 


El 


2 


J0613-02 


1.6 


0.6 


0.15 


2004.5-2009.5 






J0737-30 


44 


0.6 


0.6 


2005-2009.5 


[8] 


0.4 


J0751+18 


2.6 


3 


14 


2005-2009.5 






J1012+53 


2.6 


3.0 


2.0 


2005-2009.5 






J1022+10 


3.2 


135 


360 


2005-2009.5 






J1024-0719 


1 


0.6 


2 


2006-2009.5 






J1300+12 


9 


18 


36 


2005-2009 


[32] 


30 


J1455-33 


3.2 


1.5 


3.0 


2005-2009.5 






J1643-12 


2.1 


0.7 


0.4 


2005-2009.5 






J1713+07 


0.2 


2 


1 


2005-2008.5 


m 


2 


J1730-2304 


1.3 


1 


245 


2005-2009.5 






J1744-11 


0.9 


0.9 


0.4 


2005-2009.5 






J1824-24 


14.0 


0.3 


0.3 


1990-2008 






J1857+0943 


0.8 


1 


3 


2005-2009.5 






B1909-37 


0.3 


0.8 


0.7 


2005-2007 






B1937+21 


0.3 


0.9 


1 


2005-2007 




3 


J2145-07 


2.1 


1 


4 


2005-2009.5 




1 



of rotation matrices between planetary frames. As expected, the obtained matrices 
are consistent with results by (Folkner et al. 1994 [15]) and (Standish 1998 [41 j). 
These results confirm the mas-level consistency of the planetary ephemerides and 
check estimations of the internal accuracy reached for the Earth orbit. 

in tableware given the rotation angles obtained between planetary ephemerides 
frames and ICRF as deduced from pulsar timing and VLBI astrometry. Among 
the 18 selected MSP, only four of them have a mas level accuracy in VLBI and 
radio timing astrometry. in table [7| we give the angles between planetary frames 
only on the four pulsars. These angles are consistent with the ones obtained with 
18 objects. It seems then that no peculiar bias is introduced by this selection 
restricted to only four sources, in table are also given the rotation angles between 
the planetary ephemerides and the ICRF. The uncertainties are here about one 
order of magnitude larger than the internal accuracy estimated previously. One 
can note a significant rotation of about 10 mas in the y-axis common for all the 
planetary ephemerides. Based on the expected accuracies of the VLBI tracking of 
spacecraft from 0.5 to 10 mas, the other rotation angles in x and z-axis are not 
significant. It is no clear if the detected rotation in y-axis is real or induced by 
a the limited number of pulsars used to test the direct link between planetary 
frames and the ICRF. An increase of the number of the MSP with a mas level 
VLBI astrometry will help to validate this conclusion. We see here the efficiency 
of using radio timing and VLBI astrometry of millisecond pulsars to estimate the 
rotation matrices between the planetary frames and the ICRF. We checked that 
the procedures give similar angles between planetary ephemerides frames as the 
one already deduced by other methods (Standish 1998 [JT], Folkner et al. 1994 
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Table 7 Angles of rotation deduced from adjustment of rotation matrices following equations 
1 and 2. The angles are given in mas and the uncertainties are the formal 1-sigma deduced 
from the least squares. In the first part of the table, the angles were obtained in using 18 MSP 
observed at NRT with an astrometric accuracy better than 10 mas. In the second part of the 
table, angles obtained in using 4 MSP (J0737-30, J1713+07, B1937+21, J2145-07) observed 
in radio timing with an astrometry better than 10 mas are given. In the third part of the 
table, ICRF to planetary frame angles are obtained by using the same 4 MSP observed with 
an astrometry better than 10 mas both in VLBI and in radio timing. |15| stands for (Folkner 
et al. 1994), |4T] for (Standish 1998) 
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|15j ) and confirmed the mas level internal accuracy of the planetary ephemerides. 
More MSP should be observed in VLBI in order to increase the sample to directly 
link planetary frames and ICRF. 
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Fig. 1 Comparisons of geocentric right ascensions, declinations and distances of Mercure and 
Venus based on INPOPlOa, DE421, INPOP08 and INPOP06 
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Fig. 2 Comparisons of geocentric right ascensions, declinations and distances of Mars based 
on INPOPlOa, DE421, INPOP08 and INPOP06 
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Fig. 4 Comparisons of geocentric right ascensions, declinations and distances of Jupiter and 
Saturn based on INPOPlOa, DE421, INPOP08 and INPOP06 
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Fig. 1 MEX (line 1), MGS (line 2) and Viking (line 3) residuals based on INPOP06 (left-hand 
side), INPOP08 (middle) and INPOPlOa (right-hand side). Light colors stand for extrapolated 
residuals. 
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Fig. 2 VEX (line 1) residuals based on INPOP06 (left-hand side), INPOP08 (middle) and 
INPOPlOa (right-hand side). Light colors stand for extrapolated residuals. 
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Fig. 3 Mercure (line 1) residuals based on INPOP06 (left-hand side), INPOP08 (middle) and 
INPOPlOa (right-hand side). Light colors stand for extrapolated residuals. 
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Fig. 4 VLBI de Venus (line 1), Mars (line 2) and Jupiter (line 3) residuals based on INPOP06 
(left-hand side), INPOP08 (middle) and INPOPlOa (right-hand side). Light colors stand for 
extrapolated residuals. 
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Fig. 5 Cassini tracking residuals based on INPOP06 (left-hand side), INPOP08 (middle) and 
INPOPlOa (right-hand side) given in right ascension (line 3), declination (line 2) and geocentric 
distances (line 1). Light colors stand for extrapolated residuals. 
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Fig. 6 Cassini VLBI residuals based on INPOP06 (left-hand side), INPOP08 (middle) and 
INPOPlOa (right-hand side) given in right ascension (line 3) and declination (line 2). Light 
colors stand for extrapolated residuals. 
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Fig. 1 Optical residuals of Jupiter in right ascension (line 1) and declination (line 2) and of 
Saturne (lines 3 and 4) based on INPOP06 (left-hand side), INPOP08 (middle) and INPOPlOa 
(right-hand side). Light colors stand for extrapolated residuals. 
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Fig. 2 Optical residuals of Uranus in right ascension (line 1) and declination (line 2) and of 
Neptune (lines 3 and 4) based on INPOP06 (left-hand side), INPOP08 (middle) and INPOPlOa 
(right-hand side). Light colors stand for extrapolated residuals. 
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Fig. 3 Pluto residuals in right ascension (lines 1 and 3) and declination (lines 2 and 4) deduced 
from optical observations (lines 1 and 2) and stellar occultations (lines 3 and 4) based on 
INPOP06 (left-hand side), INPOP08 (middle) and INPOPlOa (right-hand side). Light colors 
stand for extrapolated residuals. 
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